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Abstract 

We discuss the orbital angular momentum of partons inside a longitudinally polarized proton in the 
recently proposed framework of spin decomposition. The quark orbital angular momentum defined by Ji 
can be decomposed into the 'canonical' and the 'potential' angular momentum parts, both of which are 
represented as the matrix element of a manifestly gauge invariant operator. 



1 Introduction 

Recent polarized beam experiments and global QCD analyses suggest that the contribution 
of the gluon helicity AG to the spin of the proton is rather small [IJ. This observation, to- 
gether with the inexorable fact that the quark helicity contribution AS is also small (less than 
30%), lead one to suspect that the key to understand the proton spin puzzle is the orbital an- 
gular momentum (0AM) of quarks and gluons. However, progress in this direction has been 
hindered by a number of difficulties in measuring, and even defining the 0AM. So far, the only 
well-recognized, gauge invariant definition of the quark 0AM is the one by Ji [2J which can 
be measured, indirectly, as the difference between a certain moment of the generalized par- 
ton distribution and AS. Although generally accepted, this approach may be criticized on the 
basis that the corresponding operator is not the 'canonical' one that satisfies the fundamen- 
tal commutation relation of the angular momentum operator in quantum mechanics. Efforts 
to improve upon this point have led Chen et at. to propose a completely new decomposition 
scheme of the QCD angular momentum tensor BI4II which has triggered a flurry of activity 
lately I15I6I7I8I9I10I1 1I12I13I14L However, the issue still remains very controversial, and the 
overarching impact of this new formalism as well as its practical usefulness in phenomenology 
are yet to be clarified. 

In this work, we investigate the quark 0AM along the line of our previous work ifTOl which 
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we view as the proper rendition of the formalism II3I4II in the context of high energy QCD. 
We shall show that one can represent the canonical 0AM as the matrix element of a mani- 
festly gauge invariant operator which turns out to be equivalent to that obtained in the Wigner 
distribution approach [fT6ll . This paves the way to measure the canonical 0AM experimentally 
or numerically on a lattice, and thus helps to mitigate the criticism that the matrix elements 
defined in Q do not have known physical measurements [fTTII . Actually, in the gluon helicity 
sector Ref. IfTOl has already shown how one can reconcile the gluon helicity defined in [3 1 with 
AG which is measurable. We now extend this finding to the 0AM sector. 



2 Decomposition of the QCD angular momentum operator 

The main idea of [3 4J is that one can achieve a complete, gauge invariant decomposition of 
the QCD angular momentum operator by identifying the 'physical' and 'pure gauge' compo- 
nents of the gauge field 

A'^ = v4phys + ^pure ' (1) 

j^fj-f ^ f)fj. Af _ f)" Ai^ J^in\A^^ A" 1— n (1\ 

-'pure ^pure "-^ ^pure '"i' Impure' ^purej ^i V^^ 

which transform differently under gauge transformations 

^phys ~^ ^^^phys^ i 

^ U^A^^,,,U - -U^d^U . (3) 

The QCD angular momentum tensor M'^'^^ can then be written as a sum of the helicity and the 
orbital angular momentum of quarks and gluons. In the 'covariant' form [i5i6J useful for high 
energy experiments, the original proposal by Chen et al. [.3.4,1 reads 

<:lspin = -^e"'^'^^757.^, (4) 

<ul.o,bit = ^Yix'^tDl^^., - x\D;^Ji^ , (5) 

gluon-spin a ^phys a ^phys ' 

-^^uon-orbit = ^a"{^''i^pure^^a^^)a — 3;^(-Dpure^?^')a) • C^) 

where -Dpure = d'^ + igA'^^^^^, and a, 6 = 1, 2, ■ ■ ■ , 8 are the color indices Q Using the transfor- 
mation rule dl]), it is easy to check that each of the above components is gauge invariant. Such a 



^ Our convention is e'^i^s _ = _^^o^i^2^3 ^j^^j^ ^j^g light-cone coordinates = 

-^(x'^ ± x^) and denote the transverse coordinates with latin indices xt = {x*}, = 1, 2). The 

two-dimensional antisymmetric tensor is defined as e'-^ = — e+^*-^ , e^^ = eu = —e^^ = 1. 
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complete decomposition goes beyond Ji's framework in which the gluonic part cannot be sepa- 
rated into the helicity and orbital parts. The price to pay, however, is that the decomposition is 
not local in the sense that A'^^^^^ is in general nonlocally related to the total Af^. Moreover, it is 
not entirely covariant, either, because Aphys actually depends on the frame as we shall soon see. 
An alternative decomposition of the orbital part was suggested by Wakamatsu O 

M',ZL,it = ^I'ix'tD' - x\D^)i^ , (8) 



_ptJ.a( v(r)\ /iphys\ _™A/t-)1^ /IphysN \ 
-"^^gluon-orbit -'a l-^pure^Q )a l-^pure^a )a) 

HDcF^^aix^AX, - J . (9) 

The second term of ^ is gauge invariant on its own. Using the equation of motion DaF^^ = 
g^'yH'^ip, one sees that it accounts for the difference between -Dpu^ in ® and in ([8]). 

There is no consensus as to which definition, dS]) or ([8]), is more appropriate for the quark 
orbital angular momentum. To some extent, it is a matter of choice. It is ([5]), but not ([8]), that 
is compatible with the (equal-time) canonical commutation relation of the angular momentum 
operator L x L = iL, 

L = X X impure 5 X = (X^, X^, X^) . (10) 

The pure gauge condition Q is crucial for this. ([5]) may thus be called the canonical angular 
momentumifj On the other hand, ([8]) is the same as Ji's definition [2] and is accessible from the 
analysis of the generalized parton distribution (GPD), whereas it has not been known how to 
measure (|5]). 

We will derive an explicit expression of the canonical angular momentum ([5]) in terms of a 
manifestly gauge invariant operator whose matrix element is, in principle, related to experimen- 
tal processes or observables in lattice QCD simulations. For this purpose, one must specify what 
Api^yj is. There are several proposals for A'^^^^ in the literature [I4I7I10I13II . These definitions are 
not equivalent as suggested by the work of Ref. JH which showed that they give different values 
of the gluon helicity (the proton matrix element of Here we employ the one proposed in 

m 

^pV(^) = - / dy-ICiy - x)W^^F+^{y-, x)W-, , (11) 
where we use the notation x = (x"*", x*) from now on. W is the Wilson line operator 

= Vexp i^-W A+{y'-,x)dy'-^ , (12) 



^ Throughout this paper, we associate the term 'canonical' with the operator iDp^re instead of the usual 
id^^. The former is actually the gauge covariant generalization of the latter without affecting the commu- 
tation relation. 
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in the fundamental representation. The superscript '— ' denotes that the path ordering is in the 
direction. /C(y^) is either ^e(y~), 0{y^) or —9(—y^), depending on the boundary condition 
at x~ = ±00 in the light-cone gauge ^4+ = oO] The pure gauge part Apme is 

where W±oo = V exp (^—ig A{±oo, x') ■ dx'^ is the Wilson line in the spatial direction at 
x" = ±00. It represents the residual gauge symmetry of the light-cone gauge = 0, and is 
fixed by specifying the boundary condition of the gauge field at x^ — t- ±00 mentioned above. 
It has been shown in [10] that (fTT)) and (fT3l) are a viable decomposition of the total gauge field 
A'^. We wish to stress that this particular choice is singled out among others by the criterion 
of measurability: The corresponding gluon helicity coincides with the usual gluon helicity AG 
that has been measured in experiments. 

Note that the definition (fTTl) already selects a particular frame — the infinite momentum frame 
where the partonic interpretation of hadrons is clearest. As emphasized in 0151 . the decomposi- 
tion of spin into the helicity and the orbital parts cannot be made entirely covariant, but depends 
on the frame of reference. 



3 Potential angular momentum 

We now focus on the orbital angular momentum of quarks inside a longitudinally polarized 
proton. It is given by the forward matrix element of the fiuX = +ij component of Q or ([8]) 

_ 1 {PS\Jdx~d'xTM^2^_,,,,\PS) 



2P+ (27r)353(o) 

(xHDL... - xHDi...^'^b\PS) 

(14) 



1 (PS| / dx-d-'xT ^^^{xHDl,^^ - xHDi,J^\PS) 



2P+ (27r)353(0) 



1 {PS\Jdx-d^XTM';Xor,-APS) 

^ ^'~2P+ (27r)353(o) 

1 {PS\\dx-d^XTi)i^{x'%D^ -xHD')^\PS) 
~2P+ (277)353(0) 



(15) 



where = —5^ = (the proton mass squared) and (27r)^5^(0) = J dx d^xr is the 
momentum space delta function. The longitudinal polarization means = (5+, 5*", St) ~ 



3 From the viewpoint of the PT (parity and time-reversal) symmetry which will be crucially used below, 
it seems that the choice /C(y^) = ^e{y~) is the most natural and convenient one, although the difference 
does not matter in the end. 
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(5*+, 0, Ot). We first observe that, since Ai^^^^ = in the light-cone gauge = HOll, Lchen 
is actually identical to the Jaffe-Manohar (JM) definition [|18|] of the quark orbital angular mo- 
mentum 

1 {PS\S dx-(fxTi}i^{x^id'^ - xHd'^)i)\PS\c 
^™ " 2P^ (2vr)353(o) • ^^^^ 

The subscript lc means that the matrix element is evaluated in the light-cone gauge. In other 
words, Lchen is the generalization of Ljm to arbitrary gauges (see, also, [[5]|). We thus unify the 
notations Lchen = -^jm = -^can by introducing the canonical orbital angular momentum Lean, 
and write 

L^i = Lean + -^pot ; (17) 

where the so-called potential angular momentum II19I5II is, with our choice of A'^^-^^, 



-x^V^(x)7+ / drlC{y- - x-)W-ygF+\y-,x)W-^i^{x)]\PS) . (18) 



Now consider the non-forward matrix element 

1 



2P+(27r)353(o) 



P'S'\ J dx^d^XTx'Hx)j+ J dy-K{y- -x^)W-ygF+\y-,x)W-Mx)\PS) 



where P^ = (P^ + P'^^) /2 and the momentum transfer will be denoted as A'^ = P'^^ — P^^. The 
explicit factor can be traded for the derivative with respect to A* such that 

e'^Vt (19) 

= 2^iTo{zM^^^'^'l^(°)^V dy-lC{y-)W^ygF^'^{y-)W-,m\PS) - {^ ^ j) 

Parity and time-reversal (PT) symmetry tells that the following parametrization of the matrix 
element is possible 

(P'^>(0)7+ / dy-JC{y-) gF^\y-)W;,m\PS) = ze'^A^S^h{0 + ■ ■ ■ , (20) 

where 5 = (5' + S")/2 and,^ = —A+/2P+ is the skewness parameter. [The dependence of /i(^) 
on the renormalization scale is implicit. We also suppress the dependence on A^ ^ —A"^ since 
it is of higher order. Similar comments apply to other distributions defined below.] This leads to 
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^pot = /^(0)|^. (21) 

(|20l ) thus defines the potential angular momentum as the matrix element of a manifestly gauge 
invariant operator. 



The quark-gluon mixed operator that appears in the matrix element (|20l) is familiar in the 
context of the twist-three mechanism of the single spin asymmetry (SSA). Let us pursue this 
analogy and consider the following non-forward matrix element 

x{P' S'm-z~ /2h+WZ^ QF^^{y~)W-Mz- /2)\PS) 



y 

-L,>^'^P-S,P„^{x^,X2,i) + -^e^'^^-5,A.$(xi,a;2,0 + • • • • (22) 



By symmetry considerations, it follows that \l/(xi,X2,0 = ^(2^273^17 —0 ^^^^ $(xi,X2,0 = 
— $(a;2, xi, —0)- In the forward limit, and in the transversely polarized case = 5^8^, only the 
\[^-term survives. The function ^'(xi, X2, 0) plays the cental role in the so-called soft gluonic 
pole mechanism of the SSA II20I21II . In the longitudinally polarized case 5^ ^ 5^5*+, the ^- 
term vanishes for the relevant component jiu = +j. By performing Fourier transformations, we 
find 



{P'S'imi^ J dy-lC{y-)WoygF+^{y-)W;,m\PS) 
= iP^ J dxidx2 /C(xi — X2)T^-^(xi, X2) 

= ie^''S+Ak J dXdxJC{x)^{X, x, 7 (23) 
where we switched to the notation X = x = xi — X2. The kernel is 

JC{x) = p.v.i = i + , (24) 



X 2 \x + ie X — ie 
in the case IC{y^) = |e(l/") and 

]Cix) = , (25) 

X ±ze 

in the cases /C(?/~) = ±9{±y^). Comparing with (l20l) . we obtain an alternative expression for 
the potential angular momentum 

Lp^t= [ dXdxlC{x)^{X,x,0). (26) 
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Note that, since 0, 0) = 0, different choices for K, lead to the same result, as they should. 



4 Canonical angular momentum 

Next we exploit the relation between the twist-three approach to the SSA and the approach 
based on the transverse momentum dependent distribution (TMD)0 In the longitudinally po- 
larized and non-forward case, we define 

/(x, qr. A) ^ / T^^^^"""""^""^" (27) 
x(P'^'|V;(-^-/2, -;^t/2)7+WZ, ^J^l^ ^W^^ ,i^{z- /2,zt/2)\PS) , 

2 '-^ 2 ' 2 2 

where is the Wilson line in the transverse direction at x~ = ±oo. In the forward case 
A = 0, the matrix element (|27l) reduces to the usual TMD. As is well-known in that context, 
there is freedom in choosing the path connecting the points (— z~/2, —zt/2) — t- (-2^/2, zt/2). 
Using the Wilson line that goes to future infinity and then comes back —z^/2 — )• +oo — )■ z^ /2, 
one takes care of the final state interaction. The TMD in this case is relevant to the semi-inclusive 
DIS (SIDIS). The other case —z~/2 — )■ — oo — )■ z~ /2 includes the initial state interaction 
relevant to the Drell-Yan process. For the present purpose, one may as well take the average of 
the two cases. 

The relation between (|27l) and (|20l) is revealed by taking the second moment of / in qt [12311^ 

F\x,A)^ j d^rqTfix^qT.A) 

-{P'S'\i,{-z-/2)^+ 1 dy- (lC{y- - z' /2) + lC{y- + z' /2)) 

xW-^QF^\y-)W-Mz- /2)\PS) \ , (28) 



2 ' 



where the kernel /C is in one-to-one correspondence with the choice of the Wilson line path in 
(l27l) . In the forward case, (l28l) vanishes for the longitudinal polarization by rotational symmetry 
in the transverse plane. In the non-forward case, however, the following structure 



^ The following discussion is similar to the works of Ref. II22I16II . We improve upon these works by 
fully taking into account the gauge field and the issue of gauge invariance. 

^ Cf. Eq. (39) of ll23l . Via partial integration, qJ^ is replaced by the spatial derivative d^. When acting 
on the Wilson line, it brings down the factor d^A^ = D~^A^ + which reduces to the two terms in 
(EHl). 
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fix, QT, A) ~ ^e^-^^S+qT^AJix, , ^, . q^) , (29) 



is allowed, so that (1281) is not necessarily zero. Note that, because of an extra minus sign from 
Aj under the PT transformation, the function / does not change signs when changing the 
directions of the Wilson line0 in contrast to the known sign flip of the spin-dependent TMDs 
in the SIDIS and Drell-Yan reactions [|24l . 
We then take the first moment in x 

J dxF\x,A) = ^l^^ip's'imi'-i^' -^'5')m\ps) 



-^\\{p's'\mi\^^' - ^'Dlm\ps) + {p's'\mi^K^y.m\ps) 



S-{P'S'\ml^{^^l.r. - i^^MmPS) . (30) 



The matrix element (1201) indeed appears in the second term of the first equality, but somewhat 
remarkably, it is absorbed by the covariant derivative in the first term. The final expression 
features precisely the 'pure gauge' part of the covariant derivative -Dpure- Differentiating with 
respect to A, we arrive at 



d 



o f 

= lim —T-— / dxF^x, A) 
A^o id A' J ^ ^ 

d r 

= j dxd'qr g| ) . (31) 

This is a formula relating the canonical 0AM of quarks to the matrix element of a well-defined, 
manifestly gauge invariant operatorQ The final expression in (|3TI) agrees with the 0AM con- 
structed by Lorce and Pasquini lfT6l from the Wigner distribution (neglecting gauge invariance). 
We have thus established a gauge invariant link between the Wigner distribution approach and 
the spin decomposition framework of Chen et al. 

Similarly, for the gluon orbital angular momentum one can define 



^ More precisely, f{x, q^, ^, At ■ qr) +f{x, q^^ —At ■ qx) after changing the directions. But 
the difference is immaterial in the limit A ^ 0. 

We note that a possible UV regularization of the integral (|3T| ) and its evolution require an entirely 



separate analysis. 
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g{x, qr, A) ^ / ^^^^e"^"— — (32) 
X {P'S'\F^-i-z-/2, -zt/2)WZ, ^^WT^ ,x>V±oo,f ^^.'^^^^2, zt/2)\PS) , 

2 '^"^ 2 ' 2 2 

where Aphys is as in (fTTI) . and now the Wilson lines are in the adjoint representation. Under the 
PT transformation, one gets the same operator back (up to the direction of the Wilson line) 
provided the skewness parameter ^ vanishes, which we assume here. In deeply virtual Compton 
scattering (DVCS), ^ = corresponds to the elastic scattering of the photon. 
Proceeding as before, one finds the double moment in x and qt 



(33) 



The matrix element of (|7]) thus becomes 

1 {PS\Jdx~d^XTM^,2n-orJPS) d r ^ 

2P^ (2^^) = WA^ J ^) 

= e''^\j dxd\Tql~g{x,ql), (34) 

where we parameterized, as A — t- 0, 

g{x, qr, A) = l^e"-'^ S^qT^Afgix, + • • • . (35) 



5 Conclusions 

We have shown that the physical part of the gauge field Aphys proposed in IfTOl leads to well- 
defined expressions for the canonical and potential angular momenta in terms of the matrix 
element of certain gauge invariant operators. If one defines AG as the gluon helicity, consis- 
tency requires that (|3T1) is the corresponding canonical angular momentum. Now that we can, 
at least in principle, measure Lji, Lean and the difference Lpot = Lji — Lean separately, it seems 
more legitimate to call Lean, or equivalently the 0AM from the Wigner distribution II16I25L 
the quarks' genuine orbital angular momentum since it satisfies the fundamental commutation 
relation. 

Regarding measurability, it should be possible to compute (|3TI) in lattice QCD simulations 
as in the case of the ordinary (forward) TMD (see, e.g., [|26l). Since there is no sign flip in (ITTI) 
when changing the directions of the Wilson lines, the matrix element may not be very sensitive 
to the choice of the path. If so, and if one is not interested in the x-dependence, one may first 
integrate over x in (1271 ) and connect the quark operators (at the same value of z~ = 0) by 
a purely spatial Wilson line in the transverse plane. This avoids the introduction of lightlike 
Wilson lines on a Euclidean lattice which appears to be a vague issue. For the gluon orbital 
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angular momentum (|34l) with g defined in (|32l ). one has to deal with lightlike Wilson lines even 
after the x-integration. 

Finally, from the experimental point of view, the matrix elements such as (|22l) are related to 
the twist-three GPDs [27 28J which are hard to extract. It would be interesting to see if there 
are processes in which these functions contribute to the cross section at leading order as in the 
single spin asymmetry. 
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